Gaps tunable by electrostatic gates in strained graphene 
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We show that when the pseudomagnetic fields created by long wavelength deformations are ap- 
propriately coupled with a scalar electric potential, a significant energy gap can emerge due to the 
formation of a Haldane state. Ramifications of this physical effect are examined through the study 
of various strain geometries commonly seen in experiments, such as strain superlattices and wrinkled 
suspended graphene. Of particular technological importance, we consider setup where this gap can 
be tunable through electrostatic gates, allowing for the design of electronic devices not realizable 
with other materials. 



I. INTRODUCTION 

Graphene^ is a material whose unique properties are 
a fascinating challenge in both fundamental and applied 
sciences. The basic properties of its electronic structure 
are chirality, electron-hole symmetry, and linear gapless 
energy spectrum, that is, charge carriers in graphene 
are massless Dirac fermion^. In addition, corrugations 
and topological defects create gauge (pseudomagnetic 
and even pseudo-gravitational) fields acting on electron 
states^. 

Recently, a novel state of matter, a quantum Hall in- 
sulator without a macroscopic magnetic field (Haldane 
stat^, has spawned the interest in unusual topological 
properties of band structures, leading to the prediction 
of topological insulators in two and three dimensionJ^KSl^ 
It was understood afterwards that such Haldane state 
can be realized in a graphene superlattice by a suit- 
able combination of scalar and vector electromagnetic 
potentials'^. A gap opens in the electronic spectrum, 
turning graphene into a quantum Hall insulator with pro- 
tected chiral edge states. Since, long waveleng th str ains 
in graphene induce a pseudomagnetic gauge fielcP^'^, the 
combination of strains and a scalar potential should too 
open a gap in graphene^^. The latter, if realised, would 
be of general interest. This is the subject of our work. 

Energy gap in graphene is crucial for many applica- 
tions. Its realization remains a challenging problem, 
as the transformation of electrons into holes i.e. Klein 
tunneling'^, is an significant obstacle. At present, all 
known ways of gap opening have a detrimental effect on 
the electron mobility. In both biased bilayer and chemi- 
cally functionalized graphene, one arrives at a disordered 
semiconductor with Mott variable range hopping mobil- 
itjli^HlII In graphene nanoribbons, the mobility is typ- 
ically several order of magnitude smaller than in bulk 
graphene'^. Our approach provides an attractive route 
to circumvent these limitations and might also allows for 
the design of new electronic devices. 

This paper is organized as follows. In Sec. II, we 
present a general argument for the gap opening due to 
combination of strains and a scalar potential. Sec. Ill 




FIG. 1: Two strained graphene lattice configurations (upper 
panel) and their corresponding pseudomagnetic field (lower 
panel) are studied in this paper i.e. (a) strain superlattice 
and (b) wrinkled graphene. 



considers specific realization of this effect in strain su- 
perlattices and wrinkles. Fig. [l] provide an illustrations 
of these two strain geometries. In particular, we consider 
physical setup where this gap can be tunable through 
electrostatic gates. Sec. IV considers related effects such 
as Fermi velocity renormalization, topological defects and 
interplay with magnetic field. 
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II. 



GENERAL ARGUMENTS 



Strain induces scalar and gauge potentials in 
graphene^. Fig. [l] illustrates the strain induced pseudo 
magnetic field for a strain superlattice and wrinkled 
graphene. In terms of the strain tensor, the scalar and 
vector potentials ard^^^^^ 

V{r) = g[ua:x{r)^Uyy{r)] 



[uxx{r) - Uyy{r 

Li 

2^Uxy{r) 



(1) 

Li 

where g ^ 4e\Sl, f3 = -dlog{t)/dlog{a) ^ "P^, t ^ 3eV 
is the hopping between tt orbitals in nearest neighbor 
carbon atoms, and a ^ 1A2 A is the distance between 
nearest neighbor atoms. 

The Hamiltonian for Dirac fermions in graphene, in- 
cluding a scalar potential and gauge fields due to strains 
is given by 

H = Ho + Ha + Hv 
1-Lo = vf {icFxTzdx + icFydy) 
Ua = -Vf [crxAx{T) + rzayAy{T)] 
Uv = V(v) 



(2) 



where is the Fermi velocity, and a^-. 



and 



are Pauli matrices which operates on the sublattice and 
valley indices. 

Using perturbation theory in both scalar and vector 
potential, we obtain a self energy correction which can 
be written as a cross-term to the Hamiltonian (linear in 
V and linear in A): 



S(w,f-r') 



Ua {col - -Hoy^ Hv 
Uv {ojX + -Ho) 



2j2 



nl) -Ha 



(3) 



where the identity matrix X is made explicit. Since the 
self energy contains term proportional to axCy = icjz^ 
a gap can be opened in both valleys. Next we make 
explicit the necessary relations between the scalar and 
gauge potential for this gap opening. 

We calculate the second order diagram for the electron 
self energy (effective potential) as shown in Fig. [2k. At 
low energies, = 0, and zero wave vector k = (corre- 
sponding to the Dirac point), the form for the induced 
gap reads 



A = -Tr <^ (Jz 



Vf 



I 



Jm(y [(ka),(A£a)' 



CX 



/ 



Jm(F ^) (kxAl - kyAl) 



(4) 




V(f ) = (b) 

— - V(f)~B(f) 

V(r)~Uxx+ii^-Y 

V(f)~u^(f) 
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FIG. 2: (a) Sketch of the diagram which describes correlations 
between the scalar potential and pseudomagnetic field, see 
Eq. [4] (b) Electronic bandstructure of graphene superlattice, 
where each supercell unit (as indicated in Fig. la) contains 
40 X 40 graphene lattice units. Out-of plane corrugations 
amplitude of /lo = 2 A is used, leading to a non-homogenoeus 
pseudomagnetic field (see Fig. la) less than lOOT. Various 
scalar potentials as indicated are considered, where V{r) = 
coB{f) and V{f) = coUz{r) leads to gap opening, co is chosen 
such that max[y(f)] =0.2V. 



and A^ are Fourier 



where [• • • ] is the commutator 

components of scalar and vector potentials, respectively. 
This equation shows that the gap is induced through the 
correlations between the scalar potential and the pseudo- 
magnetic (synthetic magnetic) field, = kxA^ — kyA^. 
We characterize these correlations by the parameter C 
such that 



or, alternatively, 
lim 

|r-r'|^c 



lim(W),- = C (5) 



{V{r)B{r')) = CS^^\r-f) (6) 



The parameter C has dimensions of energy. It is roughly 
given by the value of the scalar potential times the num- 
ber of fiux quanta due to the synthetic field over the 
region where the field and the scalar potential are cor- 
related. For C 7^ 0, the integral in Eq. Q diverges as 
k ^ 0. Then, the lower limit of the integral should be 
kmin ^ A/'^F, turning Eq. Q into a self consistent equa- 
tion for A. 

In the diffusive regime, where electrons with momen- 
tum k have an elastic scattering time r^, the divergence 
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FIG. 3: (a) Finite size strain super lattice of dimensions 200nm x lOOnm (supercell as indicated) and its out-of-plane position, 
(b) The corresponding pseudomagnetic field generated, (c) Spatially resolved current density at Fermi energy e/ = and zero 
temperature, under infinitesimally small source-to-drain bias. Scalar potential is defined as V{r) = Cf — en, where cd is Dirac 
energy. Here, we assume V{r) = 0.15V. (d) Same as previous, except now the scalar potential varies linearly along the width 
direction, i.e. V{y) oc and V — ±0.21/ along the two edges. Conducting channel 2 is "switched off" with this type of scalar 
potential. 



in the integral in Eq. Q has to be cutoff at a momen- 
tum kmin such that vpkrain ^ Tj7^ . For a periodic 
superlattice, the integral in Eq. fil) has to be replaced by 
a sum over reciprocal lattice vectors, G . For graphene in 
the diffusive regime, resonant scatterer d^^ * ^^ * or substrate 
charge d^ give rise to a dependence rZ^ ex ni|k|~^, 
where rii is the concentration of scatterers, so that the 
lower cutoff in Eq. Q is k^in cx Using this cutoff, 

we can write 



where a ~ 1.4 A is the distance between nearest carbon 
atoms, which sets the scale of the high momentum limit 
in Eq. Q. Next, we extend this general argument to 
specific examples. 

III. POSSIBLE PHYSICAL REALIZATIONS 

A. Global gap in strain superlattices 

We show results for the gap opening due to correlation 
between the synthetic magnetic field and a scalar po- 
tential in strain graphene superlattice. We assume that 
the graphene layer is corrugated, as depicted in Fig. [l^. 
Depending on the lattice mismatch with the underlying 
substrate, the supercell size could vary from 10 x 10 
(e.g. iridium) to ^ 50 x 50 (e.g. boron nitride) times 
that of graphene unit cell. Here, we assumed supercell 
to be 40 X 40. The in-plane displacements were relaxed 



in order to minimize the elastic energy. Details of the 
implementation based on continuum elasticity model are 
given in Appendix A. 

After minimizing the elastic energy, the resultant 
strains lead to an underlying non- homogeneous pseudo- 
magnetic field as depicted in Fig. [l^. Fig. shows 
the calculated electronic bandstructure of the superlat- 
tice, after considering various types of scalar potentials. 
Strain alone, V(r) = 0, does not produce a gap. Al- 
though strains can induced a scalar potential of type 
V(r) cx Uxx + Uyy., the correlation between this potential 
and the synthetic magnetic field due to the same strains 
is also zero. On the other hand, a scalar potential propor- 
tional to the height corrugation, i.e. V{y) oc leads to 
the appearance of a gap, albeit a small one. Such a scalar 
potential could be induced by the substrate through an 
existence of electric field perpendicular to the graphene 
layeil^. Graphene superlattices induced by commensu- 
ration effects between the mismatch in the lattice con- 
stants of graphene and the substrat^^^M^ will, in general 
lead to the effect considered here. The existence of a 
gap at the Dirac energy in strained graphene superlat- 
tice is also consistent with observations which show gaps 

in very c lean sam ples which are commensurate with the 
substrat^inOMa. 

In essence, strains generally induce both scalar and 
vector potentials. However, cross correlations between 
the scalar potential and the synthetic magnetic field, as 
in Eq. ([6|, vanish in many cases. In particular, the gap 
should be zero if the system remains symmetric with re- 
spect to inversion^^. A significantly large gap is obtained 
when the scalar potential and pseudomagnetic field are 
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FIG. 4: (a) Schematic of the elongated and wrinkled graphene device, of dimensions lOOnm x lOOnm. The transport is along 
X direction. In this calculation, the quasi- ID ripple along y is assumed to have a wavelength of Xw = 20nm. The gates are 
aligned underneath the graphene, with alternating potentials of diVb and spaced at distance dw = lOnm (same period as 
the underlying pseudomagnetic field), (b) The gates results in an approximately sinusoidally varying scalar potential in the 
graphene, with amplitude Vs. Energy gaps as a function of Vs is shown, for different amount of phase shift between the gates 
and underlying pseudomagnetic field, = 0, J, (c) and (d) shows the calculated local-density-of-states as function of y and 
energy along an x-cut in the middle of the device. 



perfectly correlated, as shown in Fig.|2]3 for the case when 
V{r) (X B{r). This observation is consistent with the 
general arguments that we presented in Sec. [llj See also 
Appendix A, where a general expression for A in strain 
superlattice is derived. 

Next, we examine a slightly different scenario, where 
electrostatic gates are used to engineer the scalar poten- 
tial so as to realise the correlation with the underlying 
strain-induced pseudo magnetic field. 

B. Local gap in strain superlattices through 
quantum transport calculations 

Here, we examine numerically the effect a local elec- 
tric scalar potential on the electronic transport proper- 
ties of strained graphene superlattices. The Hamiltonian 
accounting for nearest neighbor interactions between Pz 
orbit als is given b}P^, 

n = ^Vialai + Y,tijalaj (8) 

i ij 

where Vi is the on-site energy due to the scalar potential 
V{t) and tij = t(l — ^{dij — a)) is the hopping energy. 
aij is the new bond length after strain. To facilitate the 
application of various numerical techniques, the problem 
is partitioned into block slices as shown in Fig. |5j The 
retarded Green's function in Qq, the device region of in- 
terest, can then be written as (se#^^^ for general the- 
ory), 

g = {cfX -Ho-^L- ^r)'^ = (9) 



where e/ is the Fermi energy, and T>l/r are defined as 
= dL^ and T^R = tqrt^ respectively. Ql/r are the 
surface Green's function, which can be obtained numeri- 
cally through an iterative scheme^^ based on the decima- 
tion technique (see e.g.^^). Various physical quantities of 
interest such as the transmission, current /charge density, 
local density-of-states can be obtained once Q is deter- 
mined. See Appendix B for a more detailed description 
of the numerics. 

We consider a finite size strain superlattice of dimen- 
sion 200nm x lOOnm, as depicted in Fig. [3^. The cor- 
responding pseudomagnetic field is shown in Fig. ^jp. 
Transport in non-homogeneous magnetic field is domi- 
nated by bulk "magnetic" states known as snake stateJ^. 
Snake states has been observed in high mobility 2D elec- 
tron gas system, through controlled engineering of mag- 
netic field via lithographic patterning of ferromagnetic or 
superconducting thin films^^ 

Fig. ^ plots the current density due to current in- 
jection from the left contact, biased at Fermi energy of 
IbOmeV. As depicted, current flows in regions where 
B{r) ^ 0, along the direction ±V5(r) x z. Unlike the 
case of a real magnetic field, these snake states are non- 
chiral, with forward and backwards going states residing 
in opposite valleys. In the absence of short range scatter- 
ers, these states are relatively protected. Fig. [sJd shows 
three conducting snake channels which forms the back- 
bone for the conduction. Applying the general principle 
described in Sec. [llj we apply a scalar potential that ap- 
proximately correlates with the pseudomagnetic field of 
the middle channel to open up a local gap. Indeed, a local 
gap is opened, impeding current flow along this channel. 
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as shown in Fig. [3]i. Such a scalar potential can be re- 
alised experimentally with electrostatic side gates. This 
effect could be exploited in current guiding deviceJ^. 



C. Opening gaps in suspended graphene via 
wrinkles 



Wrinkles are common feature in very clean suspended 
graphene samples, leading to finite strains. Partial con- 
trol of these wrinkles can be achieved by adjusting the 
temperature, as in some cases, they are induced by 
the mismatch in thermal expansion coefficients between 
graphene and the substrate^^^ . We assume that the de- 
formation is described by the profile proposed iiP^, as 
illustrated in Fig. flji. The resulting synthetic field is dis- 
cussed elsewhere^'^. Fig. illustrates the accompanied 
pseudo-magnetic field. 

A sinusoidal-like scalar potential within graphene is in- 
duced by gates shown in Fig. |4^, which are tailored to 
correlate with the synthetic magnetic field induced by 
the strains. The results in Fig. ^ show that a gap is 
generated, whose magnitude is proportional to the scalar 
potential. We can estimate the gap induced by a combi- 
nation of strains which are changed by Su over an area 
of spatial scale £ and a scalar potential of value SV on a 
region of the same size. The synthetic magnetic field is 
of order B - {f36u)/{a£). Then 



C ^ p5u5V- 
a 



(10) 



Eq. (10) is consistent with numerical results obtained in 
Fig. [4]a! Fig. [4J3 and d clearly show the global nature of 
the gap generated. Effectively, the gate controlled gap 
allows the device in Fig. [4] to be operated as a graphene 
transistor. 

Next, we consider situation where correlation is less 



than perfect. From Eq. (10), we note that even for small 
variations in the strain, <C 1, the gap can be of or- 
der of the potential ffuctuations, 5V ^ if the correlations 
between the scalar potential and the synthetic field are 
maintained over long distances, £ ^ a. Fig. considers 
the case where there is a phase shift between the pseu- 
domagnetic field and scalar potential. This corresponds 
to a decreasing £ in Eq. (10). Indeed the gap reduces 



as expected. In general, the presence of a gap is robust 
against reasonable degree of local disorder, since inver- 
sion symmetry is still absent in most part, see also Eq. (4) 
and related discussions. 



A. Renormalization of the Fermi velocity 

As evident from Eq. ([3|, there are self energy correc- 
tions due to quadratic terms in the scalar and vector po- 
tentials. These terms lead to logarithmic corrections in 
the Fermi velocity via the renormalization of the residue 
of the Green function: 



as 



dE 



E=0 



In A 



c) (11) 



where A k is a momentum cutoff. 

This renormalization has been previously found in non- 
linear sigma model^^^^. These logarithmic corrections 
also inffuence the kinetic equation which describes trans- 
port processes^^, in these terms they describe a pseudo- 
Kondo effect due to interband scattering (in the Dirac 
point, the energies of electron and hole states coincide 
which provides a necessary degeneracy). Due to these 
corrections, the Fermi velocity decreases in the presence 
of scalar and gauge disorder. When this effect is studied 
simultaneously with the increase induced by the Coulomb 
interaction non trivial new phases can arise^^H^. 



B. Gauge field due to topological defects 

The sublattice and valley symmetries of graphene al- 
low for the definition of a second gauge field, which hy- 
bridizes states from different valleys, and does not com- 
mute with the intravalley gauge field due to long wave- 
length strains*^^. This field can be induced by topolog- 
ical defects, such as heptagons and pentagons. These 
defects are present at dislocations and grain boundaries, 
and they can be ordered periodically forming superlat- 
tices. If the synthetic magnetic field associated with this 
field is correlated with a scalar potential, a gap inducing 
term is generated. 

The gauge field due to topological defects has the form 



-vf 



(12) 



By using perturbation theory, as in Eq. ([3|, we obtain a 
self energy which is proportional to the cross correlations 
between the scalar potential and the gauge field, multi- 
plied by the operator TyCix. Modulated Zeeman couplings 
can also lead to synthetic fields which act on the spin, al- 
lowing for the possibility of spin gaps as well. 



C. Interplay with magnetic field 



IV. DISCUSSIONS 

In this section we discuss several issues related with 
the previous consideration and the ways of further devel- 
opment. 



The gap studied here is defined in the whole sample, al- 
though its value should be roughly inversely proportional 
to the ratio between the total area and the area where 
the synthetic magnetic field and scalar potential are cor- 
related. The sign of the gap is determined by the scalar 
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potential. Localized states will be formed at boundaries 
between regions where the gaps have different sign s, sim- 
ilar to the edge states in topological insulators^. 

A periodic magnetic field, when correlated with a 
scalar potential leads to a gap whose signs are opposite 
in the two valleys^^. A combination of this gap and the 
gap due to strains leads to gaps of different values in 
the two valleys, allowing for the control of the valley and 
sublattice degrees of freedom. For example, combined 
strain and synthetic magnetic field could be useful for 
valleytronics^^ . The realization of other synthetic fields 
might open new functionalities for graphene that cannot 
be achieved with other materials. 



V. CONCLUSION 

We have discussed a novel way in which a combination 
of long wavelength strains and a long wavelength corre- 
lated potential can lead to a gap in the electronic spec- 
trum of graphene. Such situation can occur naturally, 
because of correlations between a periodic substrate and 
graphene, or it can be engineered in a controlled way 
using electrostatic gates. Since the effect is induced by 
long wavelength, smooth perturbations, a gap can be in- 
duced without increasing the amount of scattering in the 
system. Finally, valley polarized edge states will be gen- 
erated, as the band structure of the modified system re- 
sembles the spectrum of a quantum Hall insulator. 



tic energy. The strains are 



A + 2/i 



(Al) 



where hij(k) are the Fourier transforms of the tensor 



hij{5i) = dihdjh 



(A2) 



In terms o f the s train tensor, the scalar and vector po- 
tentials areSH^ 



Using Eq. Q of the main text, we obtain 



(A3) 



/ 



a (A + 2/i)2 

^2 

klhyy{k) + k'^hxxO^) - 2kxkyhy{k) cos{36^) 



X / d'k 



(A4) 



This expression is zero, as = % + 7r. While the scalar 
and gauge potentials are correlated, their correlation does 
not contribute to the formation of a global gap. 
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Appendix A: Expression for energy gap in 
superlattice 

Strains induce scalar and gauge potentialJ^. We 
study the correlations between these potentials when the 
strains are induced by modulations in the vertical dis- 
placement of the layer, h{r). We assume that the in 
plane displacements relax in order to minimize the elas- 



Appendix B: Quantum transport methods. 

The Hamiltonian accounting for nearest neighbor in- 
teractions between Pz orbitals is given b}^, 

i ij 

To facilitate the application of various numerical tech- 
niques, the problem is partitioned into block slices as 
shown in Fig. [5] The retarded Green's function in ftp, the 
device region of interest, can then be written as (see^^^^ 
for general theory), 

g = {cfl -no-^L- ^r)'^ = A-^ (B2) 

where is the Fermi energy, and T^l/r are defined as 
= r^gLT and T^r = rgur^ respectively, ql/r are the 
surface Green's function, which can be obtained numeri- 
cally through an iterative scheme^^ based on the decima- 
tion technique (see e.g.!^. It is also useful to define the 
quantity, broadening function, T^/ji = i{^L/R ~ ^L/i?)- 
Physical quantities of interest such as the transmission 
T is given by, 

T = Tr{[Ti^]l[g]i[rnt[g^r,) (B3) 



7 



Left lead 



5 = 7 



s=j s = 



Central device domain i2- 







Right lead 



FIG. 5: Numerical approach: The graphene ribbon is par- 
titioned into block slices along the x-direction (transport) as 
indicated. Lattice interactions within each block is described 
by a. Nearest neighbor blocks interactions are represented 
by r. Device domain Qq will include the strains and scalar 
potential V{r). Left/right leads regions {O^l/r) are assumed 
unstrained and electrically doped, due to charge transfer from 
contacts. 



Energy gaps, as seen in Fig. 4b of main manuscript, is 
estimated by the onset of increase in T. The electron 
density n(r) at slice j is obtained from the diagonals 



elements of given by. 



til 



(B4) 



Local density-of-states (as seen i n F ig. 4c-d of main 
manuscript) is obtained from Eq. |B4| by simply setting 
/l = /i? = 1- Current density j(r) (as seen in Fig. 3c-d 
of main manuscript), flowing from slice j to j + 1 is given 
by the diagonal of J, given by. 



where. 



+1 



[^]r'[^?"]j+i) (B5) 



(B6) 



As apparent from Eq. (B3)-(B6), it is not neccessary to 



obtain the full matrix Q. Through commonly used re- 
cursive formula of the Green's function derived from the 
Dyson equation and the decimation technique, one could 
obtain these block elements of the Green's function, [G]]^ 
in a computationally /memory efficient manner. Details 
of this numerical recipe are described elsewher^^. 
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